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Abstract 

Cubic helimagnet FeGe has emerged as a class of skyrmion materials near room temperature 
that may impact future information technology. Experimentally identifying the detailed properties 
of skyrmion materials enables their practical application acceleratedly. Here we study the mag- 
netic entropy change (MEC) of single crystalline FeGe in its precursor region and clarify its close 
relation to the critical exponents of a second-order phase transition in this area. The maximum 
MEC is found to be 2.86 J/kg.K for 7.0 T magnetic field change smaller than that of common mag- 
netocaloric materials indicating the multiplicity and complexity of the magnetic structure phases 
in the precursor region. This result also implies that the competition among the multimagnetic 
phases can partly counteract the magnetic field driven force and establishes a stable balance. Based 
on the obtained MEC and the critical exponents, the exact Curie temperature of single crystalline 
FeGe under zero magnetic field is confirmed to be 279.1 K, higher than previously reported 278.2 


K. This finding pave the way for reconstruction of FeGe phase diagram in the precursor region. 
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I. INTRODUCTION 


Novel magnetic materials exhibit a plethora of fascinating features and aspects rendering 
them relevant for basic research as well as for potential applications. Recently, magnetic 
skyrmions, a topologically nontrivial vortex-like magnetic structure, are promising for build- 
ing next-generation magnetic memories and spintronic devices due to their stability, small 
size and the extremely low currents needed to move them[l, 2]. In particular, skyrmion-based 
racetrack memory is attractive for information technology, where skyrmions are used to store 
information as data bits instead of traditional domain walls. Most of magnetic skyrmions 
have been found in chiral-lattice magnets, e.g., MnSi with B20-type crystal structure[3]. 
Due to the lack of inversion symmetry, Dzyaloshinskii-Moriya (DM) interaction shows up in 
competition with the ferromagnetic interaction, resulting in the helical (screw) spin struc- 
ture at zero magnetic field as well as in the conical and skyrmion lattice structures in 
moderate magnetic fields. Besides MnSi, observations of such chiral skyrmions have been 
also reported in noncentrosymmetric cubic ferromagnets (Fe,Co)Si and Cu3OSeO; and in 
monolayers of Fe with a strong surface induced DM coupling[4-7]. Moreover, the skyrmions 
crystal has been also recently reported in centrosymmetric cubic perovskite manganites 
Laa-3,8r1.3,Mn3O7|8]. However, skyrmions state were mostly observed at low tempera- 
tures and under the impact of external magnetic fields. The topological spin textures of 
skyrmions in Feg5Coo,58i and MnSi show up only in a low-temperature region of T«40 
K while skyrmions crystal in Lag_2,Sr1_2, Mn3O; is observed at T=20 K[8, 9]. Therefore, 
exploration of the material possessing the room-temperature skyrmion crystal is definitely 
desirable for practical applications. 

In 2011, skyrmions crystal near room temperature was first reported in FeGe film at 
260 K by Yu et al.[10]. Up to present, there are no other reported materials having such 
high temperature. Therefore, FeGe becomes the most promising candidate for realizing 
future application of skyrmion states in spintronics. The cubic FeGe belongs to the space 
group P243 , in which the non-centrosymmetric cell results in a weak DM interaction [11]. 
For FeGe, its skyrmion phase emerges in a narrow temperature range just below 278.2 
K in the filed range from 0.15 to 0.4 kOe. Besides skyrmion phase, however, more and 
more studies show that the most enigmatic and puzzling issue of FeGe is the existence of 


precursor region which comprises the so-called A-phase and a sequence of complex phase 


transitions and crossovers in the vicinity of the onset of long-range magnetic order at Tc[12]. 
Many measurements including magnetic ac-susceptibility[13], specific heat[12], electron spin 
resonance|14], small-angle neutron scattering[15] have been performed on FeGe to investigate 
this precursor region. Nevertheless, the origin and composition of precursor region is not 
well understood and maintain a long-standing and intriguing problem in helimagnets. 
Although the precursor region only exists below Tc within a narrow temperature region 
(about 5 K), it relates with the following complex phase transition and the formation of 
long-range magnetic order at much lower temperature. Recently, the critical behavior of 
paramagnetic-ferromagnetic(PM-FM) phase transition of FeGe in the vicinity of Tc has 
been studied[16]. The results show that the macroscopical ferromagnetism is composed of 
some short-range magnetic interactions and the anisotropy of magnetic exchange exists in 
the precursor region. Therefore, the accurate determination of phase transition temperature 
(Tc) will help us to further understand the composition of precursor region and reconstruct 
the more clear phase diagram. Considering the sharp PM-FM transition near Tc in pre- 
cursor region, in this paper, our main objective is to explore the accurate phase transition 
temperature under zero magnetic field through analyzing the relationship between the crit- 
ical exponents and magnetic entropy change(MEC), which allows considerable insight into 
the beginning part of precursor region and the nature of a multiplicity of magnetic phases 
within this area. The present studies show the maximum value of MEC is 2.86 J/kg.K as the 
magnetic field reaches j49H—7.0 T, implying that the coexistence and competition between 
exchange anisotropy and magneto-crystalline anisotropy can form a stable balance which 
is not apt to change with applying magnetic field. Moreover, this result also identified the 
earlier findings that the strong competition among the multi-magnetic interactions could 
soften the magnetization in the precursor region although FeGe has a sharp PM-FM phase 
transition at Tc. Meanwhile, a new method that depends on the critical exponents and the 
MEC under the different magnetic fields was applied to confirm the centra position of the 
precursor region. The exact Curie temperature is 279.1 K rather than 282 K obtained from 
temperature dependent on magnetic susceptibility of the present single crystalline and 278.2 


K in the previous report [13]. 


II. EXPERIMENT 


A single crystalline B20-type FeGe was synthesized with high-pressure apparatus.[17] 
The magnetization was measured using a Quantum Design vibrating sample magnetometer 
(SQUID-VSM) and the no-overshoot mode was applied to ensure a precise magnetic field. 
The applied field was relaxed for 2 min before the data collection. To make sure that each 
curve was initially magnetized, the isothermal magnetization measurement was performed 
after the sample was heated well above Tc for 10 min and then cooled under zero field to 


the target temperatures. 


III. RESULTS AND DISCUSSION 


Figure 1(a) shows the temperature dependence of magnetic susceptibility (x-T) of single 
crystalline FeGe, which is measured on the warming process after zero field cooling and field 
cooling under magnetic field of 100 Oe. A sharp PM-FM phase transition occurs around 
Tc=282 K (The Curie temperatures(Tc) can be determined from the peak of dx/dT (inset 
of Fig. 1(a))) The right axis of Fig. 1(a) presents the reciprocal magnetic susceptibility 
vs. temperatures and the solid line is the fitting results by using the Curie-Weiss Law 
x=C/(T-Tew) (C is the Curie constant and Tow is the Curie-Weiss temperature). Clearly, 
the relationship between x^! and T preferably follows the Curie-Weiss Law and the fitting 
Tow =283 K is also basically agreement with the Curie temperature 282 K. Meanwhile, the 
effective PM moment (Pef —2.82V/C —2.5331g) can be also calculated from the fitting Curie 
constant C. 

Figure 1(b) shows the isothermal M-H curve measured at 4.0 K. With the increase of the 
applied magnetic field, the magnetization increases sharply, and then tends to saturation as 
fo H=0.5 T. The saturated magnetization (Ms=0.961j.,) can be obtained from an extrapo- 
lation of the high field M-H curve to jj4jH—0. According to Rhodes- Wohlfarth criterion|20], 
the degree of itinerancy can be determined from the ratio of Per, to Ms. The ratio is close 
to one for the localized moment whereas it is larger than one for the itinerant moment. 
Here, the ratio of 2.63 implies that the electrons possess an itinerant characteristic in FeGe 
sample. As we know, the interaction between spin, orbital, charge, and lattice of freedom 


would be suppressed in a ferromagnet, hence the properties of magnetic phase transition can 


be clarified with macroscopic magnetization. The inset of Fig. 1(b) presents the isothermal 
magnetization vs magnetic field around 282 K. Depended on it, the properties of phase tran- 
sition of precursor region can be uncovered. Recently, the critical behavior of bulk FeGe has 
been investigated through the isothermal M-H curves and the obtained results indicate that 
the magnetic interaction of FeGe belongs to a short-range magnetic coupling in precursor 
region[16]. 

From the Maxwell relations based on the thermodynamical theory, the magnetic entropy 
change, induced by changing the magnetic field from 0 to gH, can be evaluated from the 
temperature using a numerical approximation as follows: ASy(wH,T) = Sy(uoH,T) — 
Su(0,T) = [pre Cae aia), where 10H is the magnetic field, T is the temperature 
and oH max is the maximum of the magnetic field. In the case of magnetization measure- 
ments at small discrete Ms field and temperature intervals, AS, can be approximated 


as: ASy( Bd) = [Sr MT, ugH)d(uoH) — SE" M(Ti, uoH)d(uoH)).. Ac- 


cording to this E iiis MEC for FeGe can be calculated by utilization of isothermal 
magnetization in Fig. 1(b) and the obtained AS is given in Fig. 2(a). The ASy shows a 
maximum value near 282 K and it decreases symmetrically on both sides. The magnitude 
of AS, increases with the rise of applying magnetic field and reaches a maximum value 
of 2.86 J/kg.K for uoH=7.0 T. Obviously, compared to those of well-known magnetocaloric 
standards, the values of AS, remains relatively small. In the earlier investigations on FeGe, 
some special anomalies around Tc such as U-like transition of AC susceptibility and A-like 
variation of specific-heat transition have been found[13, 18]. Different from the general 
magnetic materials, cubic helimagnet FeGe does not promptly form a long-range ordering 
magnetic structures at T «Tc. On the contrary, at the range of ~5 K below To, there are at 
least 2-3 segmented magnetic phase embedded in the conical phase[12]. The experimental 
results reveal that these magnetic structure phases are not only temperature-dependent but 
also magnetic field-dependent[13]. T'herefore, the phase diagram in this temperature region 
is very complex. Moreover, at the range of ~5 K above To, the sample does not also show 
a pure PM state. Recent thermodynamic investigations show a remarkable anomaly of a 
sharp peak occurred at the temperature slightly above Tc[18], implying that the magnetic 
component in this temperature range may be a intermediate state instead of a single PM 
or a helical phase. Therefore, in precursor region, the applied magnetic field being a sole 


driving force is impossible to make the coexistence of a multi-magnetic phases to become an 


ordering magnetic structure. According to the classical Boltzmann relation S=kglnQ, here, 
because the systemic disordering degree 2 does not significantly reduce under external mag- 
netic field, the MEC is naturally unable to achieve a large variation. It agrees reasonably 
well with the small value of MEC obtained in present FeGe. 

Nevertheless, note that Fig. 2(a) shows a broad variation of AS, with temperatures 
around To. Such a variation covering a broad temperature range maybe beneficial for 
magnetic cooling application in some special fields. One of important parameters used to 
estimate the magnetic cooling application is the change of specific heat (AC,). The specific 
heat can be calculated from the magnetic contribution to the entropy change induced in 
FeGe by the following expression: AC,(T, ugH) = C (T, oH) — C,(T,0) = TAST Moy 
Fig. 2(b) shows the AC, vs temperature under different variations of the applied magnetic 
field calculated from the AS), data (Fig. 2(a)). One can see that the AC, changes suddenly 
around To from positive to negative value with a positive value above Tc but a negative 
value below Tc. The sum of the two parts is the magnetic contribution to the total specific 
heat which affects the cooling or heating power of the magnetic refrigerator. Beside, the 
maximum/minimum values of AC, exhibit an increasing trend with the applied magnetic 
field, and are observed at the temperatures 280/290 K in 7.0 T. 

Generally, the order of magnetic phase transition can be distinguished by the criterion 
suggested by Banerjee[19]. Negative and positive slope of H/M vs M? curves respectively 
correspond to a first-order and second-order transition. However, Franco et al. recently pro- 
posed a phenomenological universal curve which could be more accurate to judge the order 
of magnetic phase transition. It is constructed by normalizing all the ASm curves against 
their respective maximum ASy"“*[20, 21], namely, AS u=AS M/ Su"? with rescaling 
the temperature axis (0) below and above Tc as defined in Eq. 1 

0— 0. = (Tpeak - T)/(Ta — Trenk): T < Ipeak (1) 
6,=(T- Tpeak)/(Tr2 - Desk); T > Tpeak 

where T,ı and T,2 are the temperatures of the two reference points that have been selected 
as those corresponding to ASyu (T,1,Tr2)=3ASy". As shown in Fig. 2(c), all the ASu 
vs. T curves under different magnetic fields collapse into a single curve further confirming a 
second order phase transition for FeGe. In precursor region, the presence of small hysteresis 


in a magnetic cycle and multipeaks in isothermal susceptibility have been found[12, 18], 


indicating that the transition between the helical, conical, and the A-phase is of first-order 
in nature. Despite all this, these mutual transitions only occur in low field region(H<0.1 
T). As the applied magnetic field is far more than 0.1 T, these first-order transitions were 
driven to a second-order. In the present study, the maximum magnetic field was applied 
up to 7.0 T. Therefore, only a second-order PM-FM phase transition was observed. Its 
properties can be characterized by some interrelated critical exponents, £ (associated with 
the spontaneous magnetization M,), y (associated with the initial magnetic susceptibility 
Xo); ô (associated with the critical magnetization isotherm at Tc). For the observed PM-FM 
phase transition of Fig. 1(a), these critical exponents (/—0.336, y=1.352, 6=5.276) have 
been obtained by numerical calculation method[16]. Here, we can further check the veracity 
of the critical exponents and disclose their relationship between magnetic entropy change 
and critical exponents. In phase transition region, AS; can be scaled by the scaled equation 
of state (497. = lar!) [22, 23] and the AS, can be rewritten into another format: 

ASw (poH, T) = (iH) g up) (2) 
Where a and A are two critical exponents (a=2-26-y and A=6+7), £ is the reduced tem- 
perature (e- Tz28). If the critical exponents are appropriately chosen, the AS, vs. T 
curves should be also rescaled into a master curve, consistent with normalizing all the AS), 
curves with two reference temperatures (T,; and T,5). Together with the obtained critical 


exponents, the AS, vs T curves under variational magnetic fields in Fig. 2(a) were re- 


—ASy (uo H,T) 
(uoH)1-9/^ 


plotted as vs ya curves. As shown in Fig. 2(d), one can find that all the 


(HoH) 
MEC collapse on a single curve for different temperature and magnetic field. This excellent 
overlap clearly indicates that the obtained critical exponents of 6 and y for FeGe sample 
are not only in agreement with the scaling hypothesis but also very reliable. 

In the following section, we will clarify the influence of critical exponents on MEC. Fig. 
3(a) presents the magnetic field dependence of isothermal entropy change at the individual 
temperatures. Obviously, in Fig. 3(a), all the entropy changes reveal a monotonic increase 
with the applied magnetic field at different temperatures. In general, it can be assumed 
that the field dependence of the magnetic entropy follows a power law expressed as AS naz 
c a(fioH)"|24, 25], where “a” is a constant and the exponent ^n" is an exponent related to 


magnetic order and depends on temperature and magnetic field as follows: 


din(|A S, |) 
din (jig) A 


For the ferromagnet undergoing a second-order phase transition, n tends to 1 and is 


n(T, HoH) = 


magnetic-field independent at temperatures T< Tc but it tends to 2 at temperatures T7» 
Tc and reaches to a minimum value equal to 2/3 at Tc|26]. As shown in Fig. 3(b), it 
can be found that the value of n is around 1.0 and reaches 1.9 far below and above To, 
respectively, basically consistent with the universal law of n-change[27, 28]. At T~282 K, 
except the value of n(u9H—1.0 T)=0.702, other values almost concentrate around 0.665, very 
close to the theoretical data of 2/3 at To. Here, we can also further check the validity of the 
obtained critical exponents of 6 and y for FeGe. In fact, at T—To, the value of n should 
strictly obey this variation established as n(Tc)-14-5— [29]. With 8=0.336 and y=1.352, 
the calculated value of n=0.61 was completely consistent with n(10H=7.0 T)=0.607 of Fig. 
3(b). 

Besides that the value of n(T, oH) can reflect the relationship between critical exponents 
and MEC, the Curie temperature can be also deduced from it. In present, the main method 
to determinate the Curie temperature is generally based on the differential data of M(T) 
curve or the linear fitting of the inverse magnetic susceptibility vs. temperatures curves. 
However, in order to increase signal to noise ratio. the actual measurement of magnetiza- 
tion always is carried out with applying a certain magnetic field. Thus, the obtained Curie 
temperature is necessarily magnetic field dependent rather than independent. Recently, a 
new method was suggested to deduce the Curie temperature that is not only more accu- 
rate but also magnetic field independent[30]. Combining Eq. 3 and the partial differential 


Gu» 


formation of Eq. 2, the value “n” can be scaled as: 


1-a 1 _Olng( grax) 
A vA ame (4) 


Hi/s/A 


n= 


Obviously, as T=Tc, w7x=0 and n=+5; on the other hand, as AS(T=Tyear)= AS maz, 
Tay = and n=+*. So, the second term of Eq. 4 vanishes at both T=T eax and T-Tc 
and the value of n(T=Tycax) is equal to n(T=Tc). Thus, the accurate Curie temperature 
Tc can be deduced from a linear interpolation of n(T, oH) curve. Here, we choose three 
representative magnetic fields (jj9H—1.0, 3.0, 6.0 T) to deduce the Curie temperature To. 


As shown in Fig. 4(a,b,c), the Curie temperature To is about 278.7, 279.2, and 279.3 


K, respectively. Fig. 4(d) shows all deduced Curie temperatures under seven different 
magnetic fields. If we adopt arithmetic mean value to determinate the Curie temperatures, 
the accurate Tc should be 279.1 K shown as the dotted line in Fig. 4(d). There is a 
discrepancy of 2.9 K comparing with Tc=282 K determined from Fig. 1(a). 

In order to identify the reliability of To=279.1 K, we can check it from the magnetic field 
dependent refrigerant capacity power (RCP) which represents the heat to be transferred 
from the cold end (Teora) to the hot end (Thot) of the refrigerator in one thermodynamic 
cycle. RCP can be calculated by integrating the (AS - T) curves over the full width at 
half maximum using the relation RCP— ee ASyudT. As shown in Fig. 5(a), the obtained 
RCP values increase with the applied magnetic field indicating that RCP is strong field 
dependent. RCP depends on magnetic field H with the rule RCP=aH°. By using it, as 
shown in inset of Fig. 5(a), the value of b=1.206 was obtained. In fact, the value of b 
is also related to the critical exponent ó( b=1+3, where 0 is associated with the critical 
magnetization isotherm at To(Mr=7.« H!/?)). So, based on the RCP results in Fig. 5(a), 
the critical exponent 0—4.854 for single crystalline FeGe can be obtained. On the other hand, 
the value of ö can be also directly deduced from the isothermal magnetization curves at To. 
Thus, under the premise of unknown exact Curie temperature, we can choose some possible 
isothermal magnetization curves which temperatures are close to exact Curie temperature 
to deduce the corresponding 6. Then, we can compare these data with 0—4.854 deduced 
from RCP. Which discrepancy is the least, which temperature is the Curie temperature. 
Here, we choose three most possible temperatures(T=278, 279, 280 K) to deduce the value 
of 6. As shown in the insets of Fig. 5(b, c, d), the values of 6 are 4.895, 4.645, and 5.136 for 
T=279, 280, and 278 K, respectively. Obviously, the value of 4.895 at 279 K is closest to 
the 6=4.854. Therefore, T=279 K can be referred to the Curie temperature instead of 282 
K previously determined from the derivative magnetic susceptibility, which also agrees well 
with the deduced Tc=279.1 K in this work. 

The earliest investigation about the Curie temperature of FeGe was first performed by 
Lundgren et al. and it was estimated to be To=280+2 K|31]. Whereafter, the Curie 
temperature was reported to be Tc=278.2 K by Wilhelm et al. from the measurement 
of AC magnetic susceptibility|13], lower than 279.1 K determined in the present sample. 
However, the temperature of Tc=279.1 K is very close to their reported 279.6 K where the 


specific heat results show a noticeable anomaly under zero magnetic field[12]. They proposed 


that the appearance of anomaly on specific heat prior to the long-ranged FM order setting 
in at Tc implied the formation of magnetic order phase above Tç due to spin fluctuations. 
For FeGe, nevertheless, the spin fluctuations are unobvious as the Fe ions possess a stable 
magnetic moment|[32]. Considering the fact that both 279.1 and 279.6 K are obtained under 
the condition of zero magnetic field, we think the abnormal magnetic order emerged might be 
not obvious above Tc. Recently, the detailed isothermal susceptibility results also indicate 
that the curve shapes at 279.5 and 279 K are almost the same whereas the curve at 278.5 K 
starts to show a double-peaks structures|18]. Therefore, the start location of the precursor 
region can be considered at 279.1 K and the phase boundary of phase diagram need to be 
adjusted accordingly. 


IV. CONCLUSION 


In summary, we have studied in detail the magnetic and MEC for Helimagnet FeGe. A 
second-order transition from PM to FM phase around Tc has been confirmed by the scaling 
MEC. A relative small MEC obtained in this sample indicates that the coexistence and 
competition between exchange anisotropy and magneto-crystalline anisotropy can form a 
stable balance in the precursor region while the applied magnetic field is unable to completely 
transform it into a single magnetic structure phase. Moreover, a more precise method for 
determining the Curie temperature was carried on FeGe sample. We found the actual Curie 
temperature of it should be 279.1 K rather than 282 K deduced directly from the derivative 


magnetic susceptibility. 
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FIG. 1: (Color online)(a)Left axes: Temperature dependence of magnetic susceptibility for FeGe. 
Right axes: Inverse susceptibility as a function of temperature and the straight line represents the 
linear fitting according to the Curie-Weiss law. Inset shows the the derivative magnetic suscepti- 
bility. (b) presents the isothermal M-H curve at 4.0 K. Inset shows the isothermal M-H around 
Te: 
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FIG. 2: (Color online)(a) Magnetic entropy change (-AS) plotted as a function of temperature(T) 
at different applied fields. (b) The change of specific heat ACp vs. T at different magnetic fields. 
(c) Normalized magnetic entropy change dependence of the rescaled temperatures. (d) The AS vs. 


T curves under different magnetic field were replotted as ET Su OS vS. To Hy 7x Curves. 
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FIG. 3: (Color online)(a) The AS vs. oH at different temperatures; (b) The value of n vs. T at 


different external magnetic fields. 
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FIG. 4: (Color online) Left axes: The AS vs. T for FeGe at external magnetic field of 1.0 T(a), 3.0 
T(b), 6.0 T(c); Right axes: The value of n vs. T at external magnetic field of 1.0 T(a), 3.0 T(b), 
6.0 T(c). (d) All the deduced Curie temperatures under different external magnetic fields and the 


dot line is the average Curie temperature. 
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FIG. 5: (Color online)(a) Refrigerant capacity power as a function of applied magnetic field. Inset 
shows the same plot in log-log scale and the solid line is the fitting results following RCP = aH. 
(b, c, d) present the isothermal M vs. H plot at 279 K, 280 K and 278 K. Inset shows the same 


plot in log-log scale and the solid lines are the linear fitting with Mr, = DH!/?, 
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